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Abstract 



The C- Algorithm introduced in |5( is designed to determine isochronous centers 
for Licnard-typc differential systems, in the general real analytic case. How- 
ever, it has a large complexity that prevents computations, even in the quartic 

CS ' polynomial case. 

The main result of this paper is an efficient algorithmic implementation of C- 
Algorithm, called ReCA (Reduced C- Algorithm) . Moreover, an adapted version 

y^Q , of it is proposed in the rational case. It is called RCA (Rational C- Algorithm) 

CO ' and is widely used in [l| and [2[ to find many new examples of isochronous 

r — | , centers for the Lienard type equation. 
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1. Introduction 

The use of symbolic computations and Computer Algebra Systems in the 
Qualitative investigations of ordinary differential equations becomes a standard 
way, see for instance [9(. One of the important problems in this field is the 
characterization and the explicit description of isochronous centers for planar 



polynomials vector fields. See for example 111 ]. 
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This is a companion paper of [l|, [2J which are devoted to the seeking out 
of isochronous centers for real planar polynomial Lienard type equation. These 
investigations are based on the C- Algorithm introduced in [5| and used in [6] 
for the cubic case. 

The C- Algorithm in its original form has a large computational cost. In the 
real analytic case as well as in the particular rational case the careful inspection 
of the formulas used leads to more efficient algorithms, called ReCA (Reduced 
C- Algorithm) and RCA (Rational C- Algorithm) respectively. The aim of this 
note is to give a detailed description of them. 

Consider the Lienard type differential equation 

as + f{x)x 2 + g(x) = (1) 

where / and g are defined in a neighborhood of £ R, or equivalently its 
associated two dimensional (planar) system 

x = y I 

. n . M2 . (2) 

JJ = -9(%) - f(x)y z ) 

In this paper, we assume that g(0) — 0, which insures that O is a critical point, 
and xg{x) > in a punctured neighborhood of € R, which insure that the 
origin is a center. Moreover, we suppose that g'(0) — 1, so that system J2]) is a 
perturbation of the linear isochronous center x = y,y = —x. 

In our knowledge equation (1) appears for the first time in M. Sabatini's 
paper 12J , where sufficient conditions for the isochronicity of the origin O are 
given for C 1 functions / and g defined in some neighborhood of 0. 

In the real analytic case, necessary and sufficient conditions for isochronicity 
are given by A. R. Chouikha in [5|, where Theorem Q] is proved. We use the 
same notations as in (al: 



F(x) := / f(s)ds, <f>(x) := / e F ^ds. (3) 

Jo Jo 

As xg(x) > for x ^ 0, the function £ is well defined by the relation 

\^f = f g(s)e 2F ^ds (4) 

and the condition x£(x) > 0, Vx ^ 0. Such function £ is real analytic in some 
neighborhood of and £'(0) = 1. 

Theorem 1 (Chouikha, 5]). Let f and g be real analytic functions in a neigh- 
borhood of 0, such that g(0) — and xg(x) > for i / 0. Then system |T]) 
has an isochronous center at O if and only if there exists an odd real analytic 
function h, called the Urabe function, satisfying the following conditions: 

m =9(x)e F ^\ (5) 



i + K&x)) 



and 



(x) = £(x) + / h(t)dt, 
Jo 



(6) 



where F(x),(p(x) and £(x) are defined by (©-flU)- 

Taking into account (Q|, it is easy to see that $5§ and (O are equivalent. 

This theorem leads to an algorithmic method, named C- Algorithm, which 
gives necessary conditions for isochronicity by computing the first coefficients of 
the power series expansion around of both sides of (|5|) . Then, the sufficiency 
is insured by establishing explicitly the odd Urabe function. More details about 
the C- Algorithm can be found in [a, ||| @, CO] ■ 

In particular, in [3,[l|, the authors investigate the practical applicability of 
the C- Algorithm to the following family of planar polynomial systems, which 
are perturbations of the linear isochronous center x = y,y 



-x: 



x =Po(x) +pi (x)y 

y = q (x) + qi (x) y + q 2 (x) y 2 



(7) 



where po, Pi, q ,qi,q 2 € R[x], po(0) = ^(0) = 0, pi(0) = 1, qo(0) = 0, q' (0) = 
— 1, qi(0) = 0. Actually, under some restrictions, such systems are reducible 
to Lienard type differential equation ([1} with rational functions / and g, and 
in this case the RCA described in Section [3] is much more efficient then the 
standard C-Algorithm. 



2. The C-Algorithm and the ReCA 

In [5[, the change of variable u — cf>(x) = £ + L h(s) ds is introduced. Let 
us denote by g(u) the function of u represented by both sides of §5$ 



■'<"> = OTo) =9(l)eF< *' 



(8) 



The derivative of u with respect to x (resp. £) can easily be expressed in terms 
of x (resp. £), 

*± = e F ^) and ^ = 1 + h(0- 
dx at, 

The C-Algorithm is based on the two sides computations of the derivatives 
of g(u) with respect to u in terms of £ and x, which consists in computing the 
following quantities: 



Po(0 



Pk(0 



i + MO 

dP k -i(0 



i 



d£ (1 + h(0) 
Qo(x) = g(x)e F ^ 



Q ki x) = *9a=M e -*M 

ax 



(9) 



where fc > 1 and evaluating them at 0, that is Pfc(0) = Qfe(O). As 



P fe (0) 



du k 



«=o, 



and Qfc(O) 



du k 



u=0, 



by analycity, (J5J is equivalent to the equalities Qfc(0) = Pfc(O), fc > 0. 

In the next theorem we describe the Reduced C- Algorithm, which takes care 
of significant multiplicative factors appearing in the formulas. The paper [2| is 
based on it. 

Theorem 2 (Reduced C- Algorithm - ReCA). Let f and g be real analytic 
functions defined in a neighborhood of 0, such that g(0) = 0, g'(0) = 1 and 
xg{x) > for x ^ 0. Then system (JTJ has an isochronous center at O if and 
only if, there exists an odd analytic function h such that for all fc > 0, 



P k (0) = Qfe(0) 



with the recursive formulas: 



Pk(0 



d£ 



(i + MO)-(2fc-i)ft-i(0 






(10) 



Qo(x) =fl>(») 

Q fc (x) = ^LzlM _ (fc _ 2 )/(ar)Q fc _ 1 (a ! ) 

dx 

Moreover, for all M > 0, the M first necessary conditions of isochronicity, 
-Pfc(O) = Qfc(0) /or 1 < fc < M can &e obtained by computing the truncated 
power series expansions around of Pk(£) an d Qk(x) up to the order M — fc. 

Proof. First, the computation of the first two derivatives of g(u) with respect 
to u in terms of both x and £ gives 



~g'(u) 



i + h (O - $ ^& 



di 



a+Mor 



-7-0 ( x ) + 5 0*0 / (x) , 

da; 



ff"<» = - 



(l+fc(OKT^) + 3^(l + ft(0-^ 



d 2 g(x) f dg (x) 



dx 2 



\ dx 



/W + sW 



4f(g) 

dx 



,-F(x) 



These formulas strongly suggest that the fc t?l derivative of <?(u) can be written 
both in terms of x and £ as follows: 



~9 {k \u) = 



PkiO 



(i + M0) 2fc+ 



r = ftW^ M)f(l) , 



with Pf. and Qk verifying the induction formulas (ITU1) . 

Then, we assume that this is so up to the order n — 1. Differentiating the 
two sides of the equality §( n— "^(ti) with respect to u in terms of x and £ gives 

S(«)(V) - Pn ^ - Q (x) e {l - n) F{x) 

where P n (£) and Q n (x) satisfy formulas (ITU|) as expected. 

As F(0) = and h(0) = 0, necessary and sufficient conditions are given by 

i 3 /c(o) = Q fc (o),fc>o. 

When our aim is to establish a fixed number M of necessary conditions, 
we restrict ourselves to the power series expansion around of Pk and Qk for 
1 < k < M up to the order M — ft, which is the minimal necessary truncation 
order. Indeed, by formulas ([ID]), to obtain -Pm(O) and Qm{0) it is sufficient to 
compute a power series expansion around of Pm(0 and Qm(x) up to order 
(i.e. constant terms) which require the power series expansion around of 
Pm-i(0 an d Qm-i(x) up to order 1, and so on. □ 

In the practical use of the described algorithm we are concerned with a 
finite number M of necessary conditions. When a candidate for an isochronous 
center is identified, we try to write down its Urabe function under a closed-form 
expression, and prove the sufficiency using again Theorem [T] 

3. The RCA 

In this section we restrict ourselves to systems ([2]) for which / and g are 
rational functions. For this particular case, we describe an easy to handle couple 
of polynomial recursive formulas which gives the fcth derivatives of each side 
of ([5]) . Those formulas apply in particular to systems when they are reducible 
to Lienard type differential equation. The paper [I] is based on it. 

We denote f(x) = N f (x)/D f (x) (resp. g(x) = N g {x)/D g (x)), where N f , N g , 
Df and D g are polynomials such that -D/(0) = 1, D g (0) = 1 and pgcd(Nf, Df) = 
l,T>gcd(N g ,D g ) = l. 

Theorem 3 (Rational C- Algorithm- RCA). There exists a positive integer Mq, 
such that for any M > Mq the following assertions are equivalents: 

1. the origin O of system ([2]) is an isochronous center; 

2. there exists a real analytic odd function h satisfying 

P k (0) = Q fc (0) 



for all < k < M , where 



(11) 



d£ ) v w/ v w di 

Qo(x) = N g (x) 

Qk{x) = Qk-iWDgfr) ((1 - k) ^^ + (2 - k) N f (x) 

~ kQk ^{x) d -^D f {x) + d 9^1 Dg{x)Df{x) 

Moreover, as we only need the values of the P k and Qk at 0, it is sufficient to 
compute the power series expansions of Pk(x) and Qk{x) at order M — k, i.e. 
to truncate the polynomials P k and Qk up to degree M — k. 

Proof. As in the proof of the previous theorem, those formulas are found by 
induction on k: 

„( k i( v ) - Pfc (£) _ Qfc(x) g-k)F(x) (12) 

9 w (1 + /i(£)) 2fe+1 D f (x) k D g (x) k + l ' l ' 

It remains to prove that there exists a finite Mq such that the Mq first 
conditions are sufficient. This comes from the Hilbcrt Basis Theorem, and more 
precisely the Ascending Chain Condition (see [7|) applied to the ascending chain 
of ideals I j = (-Pfc(O) - Qk(0), < k < j). Then there exists an M > such 
that Im = Imo+i = ■ ■ ■ = loo- n 



4. Efficiency of the RCA 

The original C- Algorithm which is based on (|9|) will be denoted by Aq. In 
this section we study the efficiency of the algorithms resulting from Theorems [2] 
and[3l that will be denoted by A-i, A3, A\ and A§: 

• A\ is a truncated C- algorithm. It is based on the formulas ©, for which 
we apply the truncation procedure using power series. 

• A2 is the algorithm based on the formulas (fTUf of Theorem [2] where the 
truncation procedure is applied, that is ReCA. 

• A 3 is the algorithm based on the formulas (fTU|) of Theorem [5] without 
truncations. 

• A4 is the algorithm based on the formulas (fTT|) of Theorem [3] where the 
truncation procedure is applied, that is RCA. 

• A5 is the algorithm based on the formulas ([fT]) of Theorem |3] without 
truncations. 



To compare the efficiency of the above 6 algorithms we will apply them to the 
quartic system 



-y + ai,ixy + a 2 ,ix 2 y + a 3 ^x 3 y 



y = x + b 2fi x 2 + h,ox 3 + b , 2 y 2 + b\. 2 xy 2 + b 2 ,2X 2 y 2 + b^ox 4 



(13) 



that is system (15) from [2j or system (3.1) from [l[. By standard reduction 
this system is reducible to Lienard type equation (|TJ) with 



m = 



fro, 2 + b\ t zx + b 2 ,2X 2 + ai4 + 2 a,2,ix + 3 az,\x 2 
1 — <X\ ix — a 2l \x 2 — a 3 1X 3 



and 



g{x) = (l — ai t ix - a 2 ,ix 2 - a 3 ,i2: 3 ) (a; + b 2 , x 2 + 63,0a; 3 + b ifi x A ) 



There are 18 unknowns, 9 for the a^- and bij and 9 for the coefficients of the 
power series expansion of h up to order 17 (remember that h is odd). Then, 
it is reasonable to compute the conditions (|TTj) at least up to order M = 19. 
Since the depth of the isochronous center still an open problem for system 
(|13[) . then investigations need a higher number of necessary conditions. In our 
comparative study we ask for the first 30 necessary conditions by each of the 
presented algorithms. 



Order of 
derivation 


C- Algorithm 

A 


A x 


ReCA 

A 2 


A 3 


RCA 

Ai 


A 5 


10 


230,8 


0,52 


0,0 


6,7 


0,0 


5,7 


15 


5920,2 


19,1 


0,6 


523,4 


0,4 


520,0 


20 




1168,6 


4,6 




2,6 




30 






168,7 




84,9 





Table 1: CPU time in seconds on Pentium 2,4 GHz with 4 Gb of memory 

The superiority of ReCA and RCA is obvious as well as the role of the 
truncation. The absence of values means that in that case the computations 
failed by lack of memory. 



5. Examples and comments 

Let us recall that the isochronicity problem for planar cubic systems (linear 
center perturbed by cubic nonlinearity) is still open. This fact is due to the 
huge number of parameters (14 parameters). In the same time, several recent 
works have proven the power of the algorithmic methods in the characteriza- 
tion of isochronous centers. For instance we quote the normal forms approach 
established by V. G. Romanovsky [ll| and used in several of his coauthored 



papers [J, Is |3[ . Particularly, this method has proven its performance in the 
study of time-reversible isochronous centers. Indeed, the paper [J] contains the 
complete set of time reversible isochronous centers of linear center perturbed by 
cubic nonlinearity. Hence the cubic isochronicity problem is still open only in 
the case of non time-reversible systems. 

Using Algorithm RCA, we succeeded to establish several new cubic and 
quartic isochronous centers [l|, |2|. Among others, we found in [l| three families 
of new non time-reversible cubic isochronous systems: 



x = — y — 2 b-zfixy + x 2 + 2 b 2fi x 3 \ 

y = x - 4 b 2fi y 2 -2xy + b 2 , x 2 + 4 b 2fi x 2 y + 2 x 3 J 

x = -y±2 V2xy + x 2 T 2 V2x 3 

y = x ± 8 V2y 2 -2xyT% V^X 2 =F 12 \[2x 2 y + 10 x 3 



1. 2 1. 3 

x = -y- - b 2fi xy + x + - b 2fi x 

y = x - b 2 . y 2 -2xy + b 2:0 x 2 + b 2:0 x 2 y + f 2 + - b 2 . z ) x 6 



(14) 
(15) 

(16) 



In |3| time-reversible isochronous centers of homogeneous quartic perturba- 
tion of the linear center are completely established. In [1|, using RCA we found a 
large list of new non time-reversible quartic isochronous centers. We also found 
the following family of systems, for which we conjecture it has an isochronous 
center at 0: 



' = - .'/ - ( - g - 2 62,2 j x 2 y + ( — + b 2a ) x 3 y + xy 

!J OJ/ j Qj 9 2 

x h 1 2 62 2 xy + b 2 2X y 

4 -r 4 -r g 2,2 y t 2,2 y lg 



(17) 



We were able to prove the isochronicity of the system in few particular cases, for 
instance for 62,2 £ { — jz, 0, jg}, by computing explicitly the Urabe function: 



h 



y / 2£v / 2£ 2 + 32(e 2 + 12) 
{62,2(0=^} ~ 2 (e + 4) (£ 2 + 16) 



ft 



W"(?)+V»(?)+ 3 H¥) 



{62,2(0=-^} 



2 f (L(f) + 4)( i ( 

where L = LambertW is the Lambert function (see [10|), 



^{62,2=0} (£) = 



, z±K!±WJ±H E g fc + 2 ^^T^ + 2 



(2 + e 2 )(V 4 + 2C 2 + 6 
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